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We present experimental and numerical studies of nonlinear magneto-optical rotation (NMOR) in 
rubidium vapor excited with resonant light tuned to the b 2 S 1 / 2 6 2 Pi/ 2 absorption line (421 nm). 

Contrary to the experiments performed to date on the strong D\ or P 2 lines, in this case, the 
spontaneous decay of the excited state 6 2 Pi/ 2 may occur via multiple intermediate states, affecting 
the dynamics, magnitude and other characteristics of NMOR. Comparing the experimental results 
with the results of modelling based on Auzinsh et al ., Phys. Rev. A 80, 1 (2009), we demonstrate 
that despite the complexity of the structure, NMOR can be adequately described with a model, 
where only a single excited-state relaxation rate is used. 

PACS numbers: 33.57.+c,32.60.+i,42.65.-k 


I. INTRODUCTION 

Nonlinear magneto-optical rotation (NMOR) is the 
light-intensity dependent rotation of polarization of lin¬ 
early polarized light during its propagation through a 
medium subject to an external magnetic field. Over the 
years, the effect has been extensively studied, both ex¬ 
perimentally and theoretically |l]. The research has been 
driven by the desire for a comprehensive understanding 
of the physical processes responsible for the rotation of 
light polarization, as well as by fundamental and prac¬ 
tical applications of the effect. For example, a detailed 
understanding of the generation, evolution, and detection 
of quantum states in atoms, manifested at the macroscale 
as NMOR, led to the development of techniques enabling 
the manipulation for these states (quantum-state engi¬ 
neering) [2 , 3] and methods of their nondestructive mea¬ 
surement (quantum nondemolition measurements) [4j|5]. 
The effect has also been used in investigations of the 
relaxation of ground-state coherences in atomic vapor 
EHS, resulting in refinement of the techniques enabling 
the generation of long-lived (>60 s) ground-state coher¬ 
ences [9]- On a more practical side, NMOR has found 
applications in atomic clocks m, optical magnetome- 
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ters mm , narrow-band optical filters 13], and laser- 
frequency locking systems M- An interesting area of 
application of NMOR is fundamental research. For ex¬ 
ample, the effect is used in the search for non-magnetic 
spin couplings USHH, and is proposed for experiments 
focusing on the detection of constituents of dark mat¬ 
ter/energy [20] . 

To date, most of the NMOR studies and applications 
utilized alkali vapors optically excited on the strong D\ or 
D 2 absorption lines (see, for example, Ref. [12] and refer¬ 
ences therein). In such systems, good agreement between 
theoretical predictions and experimental observations has 
been demonstrated mm- The wide range of successful 
applications of the method resulted in a demand for its 
extension to different optical transitions. In this paper, 
we study NMOR under such different physical conditions, 
i.e., we explore the effect with rubidium atoms excited to 
a higher-energy state (the 6 2 Pi/ 2 state). This results in 
more complex repopulation of the ground-state levels; in 
addition to the direct repopulation of the ground-state 
sublevels, the repopulation my occur via several interme- 
diate states (6 2 S'i/ 2 , 4 2 D 3/2 , 5 2 Pi/ 2 , and 5 2 P 3/2 ) (Fig. Jl]) 
[22] . This enables the analysis of the role of these dif¬ 
ferent relaxation channels in NMOR, while preserving 
other important characteristics, such as ground-state re¬ 
laxation time and atomic density. Moreover, the smaller 
splitting of the hyperfine levels of the 6 2 Pi/ 2 states, rel¬ 
ative to the 5 2 Pi/ 2 and 5 2 P3/ 2 states, offers the possibil- 
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ity to investigate the role of state splitting on the effi¬ 
ciency of generation and probing of ground-state coher- 


A specific goal of this paper is a comparison of the 
experimental results of the, so-called, blue NMOR, where 
excitation and probing at the 5 2 Ai/ 2 6 2 Pi/ 2 absorption 

line is performed using 421 nm light, with the results of 
theoretical calculations based on the model developed in 
Refs. [23l Our aim is to verify if the real (complex 
but closed) system can be adequately described with a 
single relaxation parameter responsible for ground-state 
repopulation as was assumed in the model. 

In addition to the understanding of NMOR in a 
more complex system, observation of the effect at the 
5 2 S'i/ 2 -s- 6 2 Pi/ 2 transition offers several interesting fea¬ 
tures. For example, due to the low strength of the 
5 2 Si/ 2 —>• 6 2 Pi/ 2 transition, the excitation only weakly 
perturbs the atomic medium. Thus, blue NMOR may 
be applied as a nondestructive probe of laser cooled 
and trapped atoms, including quantum degenerate gases. 
While it is possible to perform similar measurements on 
stronger transitions (e.g., the D\ or P 2 lines) using weak 
or detuned light, application of blue light facilitates sep¬ 
aration of the light from the cooling/trapping (near in¬ 
frared) beams, which increases precision of the state de¬ 
tection (optical readout is not affected by scattered light). 
In another application, low absorption of blue light in 
water provides an opportunity to use this transition in 
remote underwater magnetometery. 

This paper is organized as follows. In the next sec¬ 
tion, we outline the theoretical approach developed in 
Ref. [23], recalling the most important results of the anal¬ 
ysis (see Appendix 0- In Sec. [m) we describe the ex¬ 
perimental apparatus used for detection of blue NMOR. 
The experimental and theoretical results are presented 
and discussed in Sec. [IV] and the conclusions are given in 
Sec. 0 Finally, the Appendix presents the derivation of 
the scaling parameters that enables precise calculations 
of blue NMOR signals and its direct comparison with 
experimental data. 


II. THEORETICAL BACKGROUND 

Nonlinear magneto-optical rotation can be considered 
as a three-stage process, consisting of pumping, evolu¬ 
tion, and probing of quantum states of atoms, constitut¬ 
ing a medium interacting with light. In the model, light 
first pumps atoms, modifying their optical properties and 
introducing medium’s anisotropy. The axis of generated 
anisotropy is defined by the light polarization. During 
the next stage, the light-modified properties evolve. This 
evolution is caused by external fields, resulting in preces¬ 
sion of the anisotropy around the magnetic field. The 
precession is accompanied by processes intrinsic to the 
system, which leads to relaxation of the atomic polariza¬ 
tion and hence medium’s optical anisotropy. In the final 
stage, modified properties of the medium affect propaga- 



FIG. 1: Relevant energy states and transitions of Rb 
atoms excited by resonant radiation at 421.6 nm. The 
solid line indicate the excitation light, while the dashed 
lines indicate channels of spontaneous emission. 


tion of light, leading to the rotation of light polarization. 
Although in real experiments these three stages typically 
occurs simultaneously, here, without loss of generality, we 
consider them independently (in Appendix [B] we eventu¬ 
ally shows a polarization-rotation formula that combines 
all the stages). 

To describe the quantum evolution of a system, one 
may utilize the density-matrix formalism. This approach 
allows determination of a quantum state of atoms inter¬ 
acting with light subject to an external magnetic field, 
i.e., the situation that we encounter in NMOR. While an 
explicit form of the Liouville equation governing quan¬ 
tum evolution of the system is provided in Appendix [Bj 
here, we only recall the most important elements of the 
derivation of the equations for polarization rotation. 

In our approach, the light-atom interaction is treated 
semi-classically and is described by the electric-dipole 
Hamiltonian. Optical excitation is accompanied by spon¬ 
taneous emission, covered in the Liouville equation by a 
spontaneous-emission operator. Both processes are re¬ 
sponsible for optical pumping of the medium and hence 
generation of its polarization. It should be stressed, how¬ 
ever, that despite the significant number of the interme¬ 
diate states present in the system, in our model we only 
use a single excited-state relaxation rate. This reduces 
the problem to the simpler “two”-level system. The ques¬ 
tion if this assumption is correct or is oversimplification 
is one of the main problems addressed with this work. 

To conveniently describe optical pumping, one may 
use the irreducible tensor basis. This approach allows 
us to introduce several simplifications into the theoret¬ 
ical description. In particular, it can be shown that in 
low light-intensity regime, absorption of a photon by an 
unpolarized atom may only induce the three lowest polar- 
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ization moments: population, orientation, and alignment 
(see discussion in Appendix [bJ. Moreover, if light polar¬ 
ization direction is aligned along the quantization axis 
the only non-zero atomic polarization generated during 
pumping stage is the alignment population distribution 
(// 20 )), while the only non-zero moment existing initially 
in the system, moment (isotropic polarization), is 

modified [251. 

During the next (evolution) stage, the only external in¬ 
teraction present in the system arises from the magnetic 
held, i.e., the interaction Hamiltonian is exclusively given 
by the Zeeman interaction, and the relaxation. The in¬ 
teraction, in its lowest order, does not change the rank of 
the polarization moment, but it can mix components of 
the same-rank moments. Thereby, the polarization mo¬ 
ments // 2 ^ can be generated in the system during the 
evolution stage. 

Rotation of light polarization occurring during the last, 
probing, stage can be calculated using the macroscopic 
polarization of the medium. This enables calculation of 
the medium’s electric susceptibility and hence refractive 
indices for two orthogonal circular polarizations of light 
responsible for light-polarization rotation. In the consid¬ 
ered case, the only polarization moments contributing to 
the effect are the p( 21 ) and p^ 2 _1 ) polarization moments 
(Appendix[B|. Since these moments are generated during 
the first two stages, this shows the role of light (genera¬ 
tion of the p( 20 ) polarization moment) and magnetic held 
in NMOR. 

It should be stressed that our model was developed 
under the assumptions of low light intensities and small 
magnetic fields. The low light intensity ensures the ab¬ 
sence of such nonlinear optical processes as alignment- 
to-orientation conversion [26], which would result in de¬ 
terioration of NMOR-signal amplitude. The assumption 
of low magnetic held ensures that the effect arises due 
to the ground-state coherences, i.e., is not not caused by 
linear magneto-optical rotation or the so-called Bennet- 
structure effect m- 

A specihc question that needs to be addressed in the 
theoretical considerations of NMOR is motion of atoms. 
As shown in Refs. ESI 121], the atomic motion should be 
treated differently when atoms are contained in a buffer- 
gas-free uncoated vapor cell, a buffer-gas-hlled cell, or a 
paraffin-coated cell. The simplest situation is in the hrst 
case, i.e., the scenario investigated in this work, when 
the alkali-atom collisions are rare, and the most impor¬ 
tant mechanism of ground-state relaxation are collisions 
with the walls. The absence of velocity-class mixing en¬ 
sures during the interrogation period ensures that de¬ 
tuning does not between pumping and probing, i.e., light 
generates and probes atomic polarization at the same hy- 
perhne transition. Thereby, the signal from each velocity 
group can be calculated independently and then summed 
over the velocity distribution with weights given by the 
distribution. This allows calculation of polarization ro¬ 
tation from the whole ensemble. 

A formula for ensemble-averaged polarization rotation 


[Eq. (B16)] was used to calculate ah theoretical signals 
compared to the experimental data throughout this pa¬ 
per. 


III. EXPERIMENTAL APPARATUS 

A schematic of the experimental apparatus is shown 
in Fig. [ 2 J An Ar + laser (Coherent Innova-400), emitting 



FIG. 2: Experimental apparatus. P is the polarizer, 
WP denotes the Wollastone prism, A/2 is the halfwave 
plate, L is the lens, DM denotes the dichroic mirror, I is 
the iris, FC stands for the fiber coupler, PD is the 
photodiode, and PM and MM stand for 
polarization-maintaining and multimode fiber, 
respectively. 

10 W of continuous-wave (CW) green (514 nm) light, 
was used to pump a Ti:Sapphire laser (Coherent 899 ring 
laser). The laser enabled generation of 550 mW of near 
infrared (IR) light (842 nm), which was then coupled 
into a single-mode optical fiber and delivered to a differ¬ 
ent part of the building, where the main experiment was 
performed. The light coming out of the fiber was focused 
on a periodically poled potassium titanyl phosphate (PP 
KTP) crystal. Nonlinear properties of the crystal en¬ 
abled conversion of the IR light (200 mW) into 0.2 mW 
of blue (421 nm) light in a single-pass configuration f28] . 
The blue-light frequency was controlled by tuning the 
Ti:Sapphire laser. This enabled tuning the blue light in 
resonance with ah transitions of the 5 2 Si /2 —>• 6 2 Pi /2 line 
of both rubidium isotopes ( 85 Rb and 87 Rb) (see Fig. [ 3 ]). 
For the measurements of NMOR signal versus the light 
intensity, the blue light frequency was stabilized by ref¬ 
erencing the Ti:Sapphire laser to an optical cavity and 
was monitored with a wavemeter (HighFinesse WS/7). 
An absorption-spectroscopy setup, employing an inde¬ 
pendent reference cell containing rubidium vapor heated 
to 70° C, was used as an additional reference for the mea¬ 
surements. 

Rubidium (natural abundance) vapor, being a 
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magneto-optically active medium, was contained in an 
evacuated cylindrical (50 mm in diameter and 100 mm 
in length) glass cell placed inside a four-layer cylindrical 
magnetic shield made of mu-metal. The shield offered a 
passive attenuation of external magnetic fields at a level 
better than 10 6 [29]. A set of additional coils mounted 
inside the shield’s innermost layer was used to compen¬ 
sate residual magnetic fields and to apply a bias magnetic 
field along the light-propagation direction. The bias field 
was generated with a 12-bit current-output card, driving 
longitudinal magnetic-field coils, enabling generation of 
up to 1-G field. The shield layers were thermally isolated 
and the innermost layer was resistively heated, along with 
the vapor cell within it, to about 90°C, corresponding to 
a vapor density of 2.4 x 10 12 atoms/cm 3 and a Doppler 
width of the transition of about 600 MHz. 

In front of the shield, the blue light was spectrally fil¬ 
tered (the IR light was directed to the wavemeter by 
a dichroic mirror) and it was linearly polarized with a 
high-quality crystal polarizer. The laser-beam diameter 
was controlled by a iris m ■ After the cell, the polariza¬ 
tion of the light was detected with a balanced polarime- 
ter, consisting of a Wollaston polarizer and two photodi¬ 
odes. The photodiodes’ differential photocurrent divided 
by twice the sum of the photocurrents provides a rotation 
angle (p of NMOR [cp = (A — J 2 )/(2 /i + 2 J 2 ) for 
The signal was stored with a computer, which was also 
used to control the card generating the current for the 
magnetic coils inside the magnetic shield. 
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FIG. 3: Energy-level diagrams of the 5 2 5i/ 2 6 2 Pi/ 2 

absorption line in 87 Rb and 85 Rb. 




Magnetic field B (mG) 

FIG. 4: Rotation of linearly polarized blue (a) and IR 
(b) light tuned to the 85 Rb F g = 3 F e transition of 
the 5 2 5i/ 2 6 2 Pi/ 2 line and the 85 Rb F g = 3 F e 

transition of the 5 2 Pi/ 2 —» 5 2 Pi/ 2 line, respectively. In 
both measurements, light intensities and beam 
diameters were the same (1.9 mW/cm 2 and 1 mm, 
respectively), while temperature of the cell was 85°C for 
blue-light measurements and 23°C for IR measurements. 
The solid red curve superimposed on the data shown in 
(a) is the theoretical curve calculated using parameters 
extracted from the experiment (see text). 


IV. RESULTS AND DISCUSSION 

Figure [4ja) shows a typical NMOR signal measured 
with blue light. Similarly, as in NMOR observed at the 
D\ line [Fig.|4jb)], the signal is centered at zero magnetic 
field and reveals a dispersive shape. For a given set of pa¬ 
rameters, the blue-NMOR signal has amplitude of 4 mrad 
and width of 80 mG (width is defined as the magnetic- 
field difference between maximum and minimum of the 
signal.) While this rotation was measured at roughly 


85° C, 5 times larger rotation was observed on the D\ line 
at 23° C, corresponding to nearly 3 orders of magnitude 
lower concentration. In both cases, the lasers were tuned 
to the same ground state (the P = 2 state) and they had 
identical intensities. This clearly shows the difference in 
strength of magnetic-field-induced optical anisotropy for 
these two absorption lines m- 

The capabilities of blue NMOR as a weakly perturbing 
probe most prominently manifest themselves when the 
widths of resonances observed at the two transitions are 
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compared; the NMOR resonance observed at the D\ line 
is nearly two times broader than the one observed at the 
5 2 5 'i /2 —>> 6 2 Pi /2 transition. This difference originates 
from power broadening, which at the D\ line dominates 
the resonance width. The absence of power broadening 
is confirmed with the measurements of the width of the 
blue-NMOR resonance versus area-average light intensity 
fFie. [5]) [32]. As shown, the width of the simal remains 



FIG. 5: Widths of blue-NMOR resonances measured as 
a function of average light intensity for two transitions 
(F g = 3 F e - black circles and F g = 2 —>> F e - red 
triangle) of the 5 2 5'i / 2 —* 6 2 Pi / 2 line. The error bars of 
the F g = 2 —>> F e points are much larger than those in 
the F g = 3 F e case, as the former signals are much 
smaller, i.e., with worse signal-to-noise ratio. The solid 
black line is the constant value fitted to the 
experimental data with the uncertainties marked with 
the shaded area. Absence of the power broadening of 
the both resonances proves the weakly perturbing 
character of the interaction. The signals were measured 
at 85° C with beam diameter of 6 mm. The shaded area 
marks the uncertainty in the determination of the 
parameter r 1 . 

constant in the entire accessible range of light intensities. 
These results confirm the applicability of our theoretical 
model developed in a low-light-intensity regime [33]. 

Next, we compare the experimental data with the re¬ 
sults of simulations. A correct description of the problem 
requires determination of several experimental parame¬ 
ters: the Larmor frequency , Doppler broadening T£>, 
Rabi frequency fi#, and ground-state relaxation rate 7 . 
While determination of the first two quantities is rela¬ 
tively straightforward based on the temperature of the 
vapor, strength of the magnetic field, and the excitation 
scheme, the estimation of the remaining two parameters 
is more challenging. 


In the considered case of an evacuated buffer-gas-free 
vapor cell and absence of the power broadening, the 
ground-state relaxation rate is exclusively determined by 
the time of flight of atoms across the light beam. In the 
naive approach, one would expect that the rate is equal 
to the mean transverse (to the light-propagation direc¬ 
tion) velocity of atoms over the mean path across the 
light beam m- However, as shown in Ref. [35], the ef¬ 
fect of transit of atoms of different velocity groups across 
the light beam is more complex. For instance, the atoms 
that spend more time within the light beam (atoms 
with smaller transverse velocities) contribute more pro¬ 
nouncedly to the NMOR signal, causing essential narrow¬ 
ing of the signal [35] . To accommodate for this effect, we 
modified the naive formula by introducing the parameter 
& 7 so that the relaxation rate takes the form 


where r is the beam radius and v tr = ^ksT/m is the 
average transverse velocity, with T being the tempera¬ 
ture, m denoting the mass of the atom, and ks being 
the Boltzmann constant. The value of k 7 = 0.42(7) for 
a 6-mm beam was determined based on the fitting of the 
Lorentz curve to the experimental data (Fig. [ 5 ]). This 
value was then kept constant for all light intensities and 
light tunings [36]. The parameter for other beam diam¬ 
eters is given in Table [I] 

Another parameter that needs to be determined for 
the calculation is the Rabi frequency. While the Rabi 
frequency can be determined from the light intensity, the 
correct modeling of NMOR requires incorporation of the 
effective Rabi frequency. This originates from the fact 
that the light beam does not have a top-hat intensity 
profile and atoms traversing the beam experience differ¬ 
ent values of the Rabi frequency across the light beam. 
Therefore, in order to correctly describe the system, we 
introduce the parameter accommodating for the ef¬ 
fect. This parameter effectively replaces a real beam in¬ 
tensity profile with a top-hat profile, providing the Rabi 
frequency in a form 

Qr = fc—(2) 

where d is the electric dipole matrix element and Eq is 
the amplitude of the electric field of light, which can be 

TABLE I: Parameters k 7 and kn calculated and 
determined based on the experimental data for different 

apertures of the iris situated in front of the cell. 


d (mm) 

fc 7 

kg al 

L,Uxp 

1 

0.17(1) 

1.02 

1.04(1) 

3 

0.23(2) 

1.15 

1.30(2) 

6 

0.42(7) 

1.60 

1.62(7) 
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calculated from the light intensity I [I = cEq/(8tt)]. The 
parameter kn can be calculated by integrating the Rabi 
frequency over the actual beam profile (see Appendix |A| 
for more details) [37] . 

Table [T] presents the parameter k^ al calculated numer¬ 
ically based on the beam profile and iris apperture and 
the parameter k^ xp arising from the fitting of the Lorentz 
curve to the experimental data. As shown, the values of 
the parameters depend on opening of the iris (and hence 
beam diameter), as the aperture determines the profile 
of the beam. In particular, with a fully open aperture 
(6 mm), atoms interact with a beam of the Gaussian 
profile, while closing the aperture to 1 mm changes the 
profile to nearly top-hat shape. Thereby, in the first case, 
the actual beam profile strongly deviates from the top- 
hat shape and the value of kn is significantly larger than 
one, while in the second case it is nearly one. 

Determining & 7 and kn enables one to simulate the 
blue-NMOR signals. Figure |4ja) presents the experimen¬ 
tal data overlaid with results of theoretical calculations. 
The simulated curve reproduces the features of the exper¬ 
imental data including the signals’ shapes, amplitudes, 
and widths. It should be stressed that besides the pa¬ 
rameter fc 7 , the agreement between experimental results 
and theoretical simulations was achieved with parameters 
directly extracted from the experiment. 

To further check the theoretical model, we studied the 
dependence of the NMOR signals on the light intensity. 
As shown in Fig. [5j in the accessible intensity range the 
width of blue-NMOR signal is independent on the inten¬ 
sity. At the same time, the nonlinear character of NMOR 
implies the dependence of the amplitude of the signal 
on the intensity. Figure [6] shows the amplitude of blue 
NMOR signal as a function of the averaged light intensity 
measured at two transitions of 85 Rb (F g = 3 F e and 
F g = 2 -> F e ) (Fig.§. The presented data sets reveal 
different dependences on the average light intensity. For 
a given beam diameter (i.e., fixed k 1 and &q), this differ¬ 
ence is well reproduced by theory (the difference stems 
from different dipole matrix elements associated with the 
transitions). 

The dependence of the NMOR signal on beam di¬ 
ameter is shown in Fig. [7] Analysis of the ampli¬ 
tude and width versus the average light intensity re¬ 
veals some interesting features of the dependences. First 
is the difference in the slopes of the amplitude depen¬ 
dence on the average light intensity for three iris aper¬ 
tures [3.7(1) mradxcm 2 /mW, 12.5(4) mradxcm 2 /mW, 
22(2) mrad x cm 2 /mW for 1-mm, 3-mm, and 6-mm open¬ 
ings, respectively]. This originates from the differ¬ 
ent ground-state relaxation rate in the three cases [see 
Fig- 0b)] and hence different values of the saturation 
parameter k s . However, to adequately reproduce the de¬ 
pendences, one also needs to take into account the differ¬ 
ent beam profiles of the beams [getting into the satura¬ 
tion parameter via the parameters kn (Table [T])] . Second, 
the width of the resonance, although independent of the 
light intensity, depends of the iris opening [Fig.[7][b)]. In- 



FIG. 6: Amplitudes of NMOR signals measured versus 
the average intensity of light tuned to different 
ground-state hyperfine states of 85 Rb. The 
experimental data are in agreement with theory. The 
signals were measured at 85° C with 6-mm beam 
diameter. The shaded areas mark the uncertainty due 
to the limited precision of fto-parameter determination. 


terestingly, the widths are not in straightforward relation 
with the diameter of the iris opening. The experimentally 
determined widths for 3-mm and 6-mm openings are the 
same within the error bars [31(1) mG and 27(4) mG, re¬ 
spectively] , the width for the smallest opening is roughly 
2 times broader [68(5) mG]. As the blue-NMOR signals 
are not power broadened, this difference is not related 
with the effect and it stems from the beam profile [38] 
but also is a consequence of the different dynamics of the 
optical pumping for distinct beam diameters [35] . 

The final step in our studies of blue NMOR was the 
investigation of the NMOR spectra. Figure [8] shows ro¬ 
tation of polarization as a function of blue-light detuning 
for the whole 5 2 Si/ 2 —>* 6 2 Pi/ 2 line of 85 Rb and 87 Rb. The 
signal was measured at a magnetic field of roughly 20 mG 
and light intensity of ^0.1 mW/cm 2 . As discussed in 
Sec. [TTJ the absence of buffer gas and paraffin-coating of 
the walls results in the absence of “communication” be¬ 
tween atoms from different velocity classes. Therefore, 
atoms from each of the velocity classes contributes inde¬ 
pendently to the signal and the theoretical signal may 
be simulated by simply summing the contributions from 
atoms of different velocity groups. We apply this ap¬ 
proach for calculating NMOR spectra. An example of 
such simulations is presented in Fig. [8](a) along with the 
experimental data. As shown, the agreement of the simu¬ 
lations with the experimental data is good. For example, 
the calculation and measurements show that for a given 
magnetic field the strongest rotation was observed at the 
F g = 3 F e transition of 85 Rb [39]. This rotation was 
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FIG. 7: Amplitude (a) and width (b) of blue NMOR 
signal measured versus average light intensity for three 
iris apertures (i.e., for different beam diameters and 
profiles). All three experimental data sets reveal good 
agreement with the theory. The signals were measured 
for the laser tuned to the F g = 3 —)> F e transition of 
85 Rb. 


roughly 2 times larger than the amplitude of the second 
largest peak observed at the F g = 2 —>► F e transition 
of 87 Rb. The rotation ratio at these two transitions is 
roughly the same as the ratio of absorption at the tran¬ 
sitions. Much weaker signals were observed at the other 
transitions (the F g = 2 —>> F e transition of 85 Rb, and 
the F g = 1 —>> F e transition of 87 Rb.) In contrast to 
the former case, however, the difference in NMOR-signal 
amplitudes is much larger than the difference in the ab¬ 
sorption at the lines. Moreover, the signs of the rotation 
at these two lines is reversed with respect to the rota¬ 
tion observed at the F g = 3 F e and F g = 2 —>> F e 
transitions of 85 Rb and 87 Rb, respectively. The oppo¬ 


85 Rb 
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FIG. 8 : (a) Measured and simulated NMOR spectra 
recorded for a magnetic field of about 8 mG and light 
intensity of 0.1 mW/cm 2 and (b) rotation spectra 
measured for different magnetic fields and light 
intensity of 0.1 mW/cm 2 . The signals were obtained for 
light frequency scanned across all transitions of the 
5 2 5 i /2 ~^ 6 2 Pi /2 line. 


site sign stems from the opposite Lande factor of the 
two ground states. These observations were confirmed 
with the spectra measured at different magnetic field for 
slightly smaller light intensity (0.1 mW/cm 2 with respect 
to 0.2 mW/cm 2 ) [Fig. gb)]. 


V. CONCLUSIONS 

We presented results on nonlinear magneto-optical ro¬ 
tation at the 5 2 5'i / 2 —» 6 2 Pi / 2 line (421 nm) of both stable 
isotopes of rubidium. In contrast to NMOR measured on 
the D i or line, excited atoms in the considered case 
may decay through several intermediate states. Gener¬ 
ally, this modifies the dynamics of repopulation pumping, 
potentially changing the characteristics of NMOR signals 
and their dependences on various experimental parame- 





























ters. With our results we demonstrated that the theoret¬ 
ical description agrees with experimental data, reproduc¬ 
ing such characteristics of NMOR signals as their ampli¬ 
tudes, widths, and spectral dependences. This agreement 
was achieved with the model where only one excited-state 
relaxation rate was used. This is an interesting result as 
it was achieved in the system where only 15% of atoms 
return directly to the ground state, while remaining 85% 
reach the state via several intermediate states. 

The agreement between the theory and NMOR exper¬ 
iment with a buffer-gas-free uncoated vapor cell implies 
investigation of the effect in systems with velocity group 
mixing, being more sensitive to hyperfine splittings of 
the states. This can be realized in a buffer-gas-filled or 
high-temperature paraffin-coated cell. The investigations 
would further test the model developed in Ref. |23j. 

Beyond the fundamental understanding of NMOR in¬ 
duced by light coupling ground states with higher excited 
states, our results demonstrate the potential of such an 
excitation scheme for magnetometry. This would be par¬ 
ticularly interesting in remote underwater magnetome¬ 
try, due to weak absorbtion of blue light in water [40] . 


where I\ is the light intensity. In orde r to substitute 
the real light-intensity profile [Eq. (|A2| with the same 


radius and power top-hat intensity profile [Eq. A3 


needs to introduce the normalization parameter kn. This 
parameter can be determined by 


kn = 



nR 2 n 


Vo 


dr' 


(A4) 


It should be stressed that the parameter only depends 
on iris diameter and is independent of the light power, 
thus we keep it constant for all light intensites of a given 
beam diameter. 

The key experimental parameter that is required to 
calculate k q is the beam radius a. From the measure¬ 
ments of the beam profile, we determined the radius as 
2 mm. This value enabled to calculate the parameters 
kn for all beam diameters (see Table [I]). 


Appendix B: Model derivation 
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Appendix A: Calculating the parameter kn 

To estimate the parameter we consider a light 

beam of a Gaussian intensity profile 

/ = J 0 e <r 2 , (Al) 

where Jo is the maximum intensity at the center of the 
beam and a is the beam radius at which intensity drops 
1/e 2 before the iris. When the beam reaches the iris, 
the iris cuts part of the beam off, modifying its beam 
diameter and intensity profile (Fig. [9|. In such a case, 
the total light power is then given by 

l* Ra /“2 

P = 27r/o / r'e~^ r dr', (A2) 

Jo 

where R a is the radius of the iris aperture. Analogically, 
the top-hat intensity profile used in our model is 

(A3) 


Dynamics of a quantum system can be described using 
the density-matrix formalism. In the considered case of 
atoms contained in a vapor cell subject to external fields 
(e.g., light and magnetic field), the time evolution of the 
density matrix p is governed by the Liouville equation 



FIG. 9: Schematic depicting the effect of an iris on the 
light-intensity profile of light used in the experiment. 
For the iris openings significantly larger than the beam 
diameter (half-width at half maximum,) the real profile 
strongly deviates from the top-hat function. At the 
same time, the small openings result in a profile close to 
the top-hat profile. 


P = IiirRl, 
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(h = 1 , c = l) 

p=-i[H,p]-^{T,p} + A + Tr(Sp). (Bl) 

Here, H = Hq + U is the Hamiltonian of the system, 
consisting of the unperturbed Hamiltonian Hq and the 
interaction operator V, T is the relaxation operator, A 
is the repopulation operator describing density-matrix- 
independent repopulation processes such as those in¬ 
duced by wall collisions, and S is the repopulation oper¬ 
ator, which describes processes depending on the density 
matrix, e.g., spontaneous emission. 

In the context of NMOR, it is useful to present quan¬ 
tum operators, including the density matrix p, in the ir¬ 
reducible tensor basis. In that case, a quantum operator 


A can be written as 

a = Y, a 7f' t %’ ( B2 ) 

Kq 

where, is the g-th component of the irreducible 

tensor T^) with rank k (k = 0,..., 2 F and q = — ft,..., ft), 
and A K ~~ f are the expansion coefficient called polariza¬ 
tion moments (PMs) or multipoles. The subscripts F 
and F' indicate the set of quantum numbers enabling 
unambiguous determination of the state with the total 
angular momentum F and F' that the operator acts on. 
The polarization moments A 1 * - 1 - of the operator A are 
related to the set of the matrix elements Apm^'m' of the 
operator A in the | Fm) basis, where m and m' are the 
magnetic quantum number, through the Wigner-Eckart 
theorem 


A { ~ A = Tr 

FF' 


( T Sf') + 1 = {FmF’ -m'\Kq)A F m,F'. 


(B3) 


Using Eq. (B3), one can expand Eq. 


in the irreducible tensor basis, which for the ground state takes the form 


pp% = i(-i) F ‘> +F 9 +K + 1 ^ p A} 

3 9 K' ,k" ,q' ,q" ,F" l 9 9 ) 


m 


F a F n 


2 FgF" I 1 F" F' g 


(«V) _ \ A K ”q”) — A K A') 


- p 


F a F" 


T —T^ // ~ // ^ 
F"F' g ^2 F"F' g 


+ 


+ A g>, + ^ 0 3 Y, ( F M\F)p { ^,(Fi\\d\\F^)(-l) F ^ +K+1 {lp 3 * * * * p 9 ,}, (B4) 

9 9 F e ,F^>F g ,F' £ l e e J 


where (F g \\d\\F e ) is the reduced dipole matrix element given by 


(Fg\\d\\F e ) = (-l)A + /+F 9 +y [Fs][Fe] j ^ Fe | j {JgM J e) 


(B5) 


with / being the nuclear spin and J being the electron 
angular momentum. For the compactness of the formula 
above, the symbol [x\ = 2x+l was introduced. (J g \\d\\J e ) 
appearing in Eq. (B5) can be calculated from 


Rle = 




3 2 J e T 




(B6) 


where R is the branching ratio for the transition J e J g , 
7 e is the excited-state relaxation rate, and c^o is the tran¬ 
sition frequency betw een ground and excited states. The 
first two terms of Eq. (Bl) are accommodated in Eq. (B4) 


by t he e xpression under the sum. The second term in 
Eq. (B4) arises due to density-matrix-independent repop- 


I 

ulation, while the last term comes from the repopulation 
due to the spontaneous emission. 

To induce the polarization rotation of light, an atom 

first needs to be polarized/pumped with light, next the 

polarization has to evolve due to the external magnetic 

field and finally this magnetic-field-modified polarization 
has to affect the polarization of light. While in real sys¬ 

tems these stages occur simultaneously, to calculate po¬ 
larization rotation, it is useful to consider them indepen¬ 
dently. 


Equation (B4) can be used to calculate the evolution of 


the density matrix during the first pumping stage, when 
atomic polarization is generated, but also during the sec- 
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ond stage when the polarization evolves due to the mag¬ 
netic field. The difference between these two stages man¬ 
ifests in the form of the interaction operator V, which, in 
the first case, is determined by the electric-dipole inter¬ 
action between light and atoms and in the second case 
by the magnetic (Zeeman) interaction. In turn, the in¬ 
teraction operator at the two stages takes the form 

ypump = _ E . d ( B7 ) 

V evo1 = -n B, (B8) 


where E is the electric field of light, B is the magnetic 
field, d is the electric dipole moment, and /i is the mag¬ 
netic dipole moment. In the description of th e evolution 
stage, we also drop the last term of Eq. (B4) as there is 
no optical pumping. 

To calculate rotation of the linear polarization of light 
traversing the medium, one may use the wave equation 


d I 2 E d 2 E _ d 2 P 
dl 2 dt 2 n dt 2 


(B9) 


where l is the distance along the light propagation direc¬ 
tion and the relation of the medium’s dipole polarization 
P with the density matrix p is given by P = 7V at Tr (pd), 
where 7V at is the atomic density. The polarization P can 
be related to the density matrix p written in the rotating 
frame via 


P = Ft E 2R e(p Fg m g ,F e m e d Fe m s ,F g m g )e tl ' k ' r Ut \ 

( B1 °) 

where dF g m g ,F e m e is the electric-dipole matrix element 
between the ground state \F g m e ) and excited state 
|F e m e ) (the sum runs over all ground- and excited 
states), pF g m g ,F e m e is the corresponding rotating-frame 
density-matrix element (optical-coherence amplitude), 
and uj is the light frequency. Based on this approach, it 
can be shown that polarization rotation of light is given 

by 


dip 

U 


A.FUJN&- 


F n F e TTl q 7TL e 


Im 


PF g m g ,F e m e ( ^F e m e ,F g m g 


(k x e) 


(BH) 


where k is the unit vector of the wave vector k and e is 
the light-polarization unit vector. 

To make use of the formulas above, the electric- 
dipole moment dF e m e ,F g m g and the density-matrix ele¬ 
ment pF g m g ,F e m e need to be presented in the irreducible 

I 


r 


tensor basis. Then, Eq. (Bll) takes the form 


dip 

— = - 27 rcj 7 V at Im 
dt 




^_ y ^ 

e • (3 • (fc x e) 


(B12) 


where the tensor (3 is given by 




/3 = E 

F g F e Kq' q' 


(- 1 ) 


Fg-\-F e ~\-ti 


^FgFa 


-e- q >e- q »(lq , lq' / \K q'+q") 


1 1 
F F 

r 9 r 9 


F 

F P 


\(F g \\d\\F e )\*p' 


2 (k q'+q") 


F a F a 


(B13) 


with e q being the sphe rical-basis unit vectors. It should 
be noted that Eq. (B12) was derived within the first-order 
perturbation theory for the optical coherence. In the 
derivation, the ground-state coherences between different 
hyperfine states were neglected. 

To explicitly relate the polarization rotation 
[Eq. ( |B12 )] with the system parameters but to present 
it in a simpler analytic form several assumptions, arising 
from the fact that a photon is a spin-1 particle, may 
be introduced. Thereby, absorption of a photon by 
an unpolarized atom may only result in generation of 
the three lowest PMs ( p q ^ with rank k < 2 ) called 
population (k = 0 ), orientation (k = 1 ), and align¬ 
ment (k = 2). Moreover, in the case of interaction of 
unpolarized atoms with light linearly polarized along 
the quantization axis, the symmetry of the problem 


I 

indicates that only the zero component of the atomic 
polarization (q = 0) may be generated. Finally, since 
there is no preferred direction but only a preferred axis, 
as a result of light-atom interaction only the polarization 
moment p( 20 ) can be generated in the system (the 
explicit formula for p( 20 ) is given in Refs. [23] [24]), while 
simultaneously the polarization moment p(°°) is being 
modified in the interaction. 

Using the same symmetry arguments as above, it can 
be shown that rotation of light polarization can be only 
caused by the PMs w ith rank f = 1 or k = 2. In fact, 
analysis of Eq. (B13) reveals that the only PM affecting 
the polarization rotation is p ~ ^. Taking this into ac- 

F 9 F 9 

count, one can calculate the rotation of the polarization 
plane of z-polarized light traversing an atomic sample 
with completely resolved hyperfine structure (ujf 9 f' 
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7ej ^F e F' e > 7e) [H] 


A) 


dip 3 
dZ = 4*' 


E(-i ) 2r ‘ 

F g F e 


iFemi 

[I] 


( (_1)2/+2J 9 J' l 1 2 

V IW] l ^ ^ 


1 1 2 

F F F 
r g r g r e 



+ R{-l) Fg Fe 


where £q = (1 /X x dX/dt) 1 is the unsaturated resonant 
absorption length. The Lorentz profile L(w FeF ) is given 


by 


L(lU f Fa ) ~ 


(7e/2) 2 


(7e/2 ) 2 


■ C J 


(B15) 


F e F„ 


with uj r FeFg being the Doppler shifted transition frequency 
of a pair of hyperfine states (uj f F = —A co FeF -f k • v, 
where Aco FeFg is the detuning of light from the optical 
transition, v denotes the velocity of atoms, interacting 
with the light), x Fg is the magnetic-lineshape parameter 
given by 


_ (7 g/ty^Fg 

XF °~ ( 7s /2)2 + fi|y 

Here, Q F = g Fg gBB is the Larmor frequency of the 
ground-state hyperfine level F g (g Fg is the Lande fac¬ 
tor of the ground state F g and /i# is the Bohr magne¬ 
ton), 7 g is the ground-state relaxation rate, and k s = 


Fg Jg 1 } XF A L ( UJ F e Fg)] 

(2J e + 1) | F Fe Fg } | 


F F 

r 9 r 9 

F F 

1 e 1 e 



(J g \\d\\Je) 2 ^)/{igle) is the reduced optical saturation pa¬ 
rameter. 


Although Eq. (B14) concerns motion of atoms by in¬ 
troducing the Doppler shift, the equation only concerned 
atoms from a single velocity group. To properly describe 
the atomic ensemble, one needs to include the velocity 
distribution of atoms in the sample. In the case of atoms 
contained in buffer-gas-free uncoated vapor cell, polar¬ 
ization of an optically pumped atom is preserved un¬ 
til a collision of the atom with a cell wall [42] • Thus, 
atomic velocity and hence light detuning do not change 
between pumping and probing stages and light generates 
and probes atomic polarization generated in a given hy¬ 
perfine state. In such a case, the signal from each velocity 
group can be calculated independently and then summed 
over the velocity distribution with weights given by the 
distribution, leading to the ensemble-averaged polariza¬ 
tion rotation ESIEU 


D *FL\ - 3 

^0 . Ft a 


dJt/ x 


£(-D 


2 Fr 


F a F e 


[Femi 

[I] 


( (-l) 27+2J s / 112 
V l F 9 p 9 F e 


1 1 2 

F F F 

1 9 1 9 1 e 


\( Je F e I Y 

J l Fg Jg 1 J 


e A FeFg /r DxF 


X 


+ i?(-l) F «- F '(2J e + l) 


1 1 2 1 f Fg Fg 2 1 f J e F e I Y 2 \ 

Fe F e F g j\F e F e l)|F 9 J g 1 / f 

(B16) 


where 

fldiy 1 = T D (2J g + l) 

° \1 di ) RN at X 2 7 e (2 J e + 1) ’ 


(B17) 


ening. Since our experiment is performed in a buffer-gas- 
free uncoated vapor cell, we use Eq. (B16) to simulate 
blue NMOR signals. 


with A being the wavelength and T F the Doppler broad- 
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